Electrical generation of spin in crystals with reduced symmetry 
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We propose a simple way of evaluating the bulk spin generation of an arbitrary crystal with a 
known band structure in the strong spin-orbit coupling limit. We show that, in the presence of an 
electric field, there exists an intrinsic torque term which gives rise to a nonzero spin generation rate. 
Using methods similar to those of recent experiments which measure spin polarization in semicon- 
ductors, this spin generation rate should be experimentally observable. The wide applicability of 
this effect is emphasized by explicit consideration of a range of examples: bulk wurtzite and strained 
zincblende (n-GaAs) lattices, as well as quantum well heterojunction systems. 
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I. INTRODUCTION 

A considerable amount of attention has been devoted 
in recent years, both experimentally and theoretically, to 
the field of semiconductor spintronics. There have been 
numerous suggestions of advantages offered by the ma- 
nipulation of the spin degree of freedom, including the 
increased functionality of spin devices, low power con- 
sumption, integration with existing technologies and the 
fact that in spin transport quantum coherence can be 
maintained on large time scales. Semiconductors can be 
used in the same spintronic devices as metals, and have 
the additional advantages brought about by the existence 
of an adjustable bandgap and by the ability to manipu- 
late the carrier density over many orders of magnitude 
by doping, gating and heterojunction formation. More- 
over, spin accumulation in a semiconductor will generate 
a much larger voltage because the density of states at the 
Fermi energy is lower than in a metal, leading to a larger 
spin splitting 1 . 

To this end it would be extremely desirable if a prac- 
tical method existed for the efficient generation of spin 
polarization inside a semiconductor, as well as for the 
transport of spins over a sizable length scale. Although 
optical spin injection-, has been known for decades, it is 
impractical for devices, since it is not sufficiently local for 
nanoelectronics. On the other hand, spin injection from 
a ferromagnetic metal into a semiconductor requires long 
spin lifetimes 3 and is impeded by the resistivity mismatch 
between the two materials^. Progress has been achieved 
with ferromagnetic semiconductors^iSiLSiS, but room 
temperature semiconductor ferromagnetismiS has not 
yet been clearly established. Recent research has focused 
on the electrical control of spins in semiconductor9ii*i2iiii, 
including theoretical work on the possibility of observing 
a spin-Hall eRecii&lii, together with various efforts to 
generate a spin currentiSiiLiS. 

Following this energetic theoretical enterprise, recent 
experimental worki2i2S has demonstrated the detection 
of a sizable spin accumulation in semiconductors to be a 
feasible task. Motivated by these findings, in this arti- 
cle we present a general theory of the intrinsic electrical 



spin generation that occurs generically in spin-orbit cou- 
pled systems that are not inversion-symmetric. We will 
concentrate on systems where the spin-orbit interactions 
are strong. The existence of the mechanism we discuss 
has been pointed out in the previous literatur o 21 ' 22 i 23 ' 24 . 
Our approach is to be contrasted with earlier theories 
which consider directly the response of a spin polariza- 
tion to an electric field, by Aronov et alM* for exam- 
ple. When calculating the spin density directly via linear 
response it is, in general, not possible to separate the 
intrinsic generation terms from the extrinsic scattering 
effects, thus physical transparency is often sacrificed. On 
the other hand, in the present paper we recognize spin 
generation as an intrinsic effect so that it may be deter- 
mined from first principles calculations. The interplay of 
the spin generation and relaxation terms, resulting in a 
finite spin polarization, emerges in the final analysis, as 
pointed out by Edelstein 22 . So far, experiment^ has at- 
tempted to measure the spin generation rate separately 
from the spin relaxation time in the weak spin-orbit cou- 
pling limit. The idea we discuss has already been applied 
to the Rashba Hamiltonian 2 i, which we also examine be- 
low. 

Spin-orbit interactions can be important in semicon- 
ductors for several reasons. The first is the fact that the 
carriers are clustered near the band extrema around high 
symmetry points where there exist degeneracies, and the 
form of these degeneracies is determined by the spin-orbit 
interaction. Due to the fact that in semiconductors the 
carriers occupy a narrow width of fc-space the spin-orbit 
interaction can therefore play a crucial role. We will show 
in this article that, because of the fact that spin is not 
conserved, there exists a term which acts as a bulk source 
of spin generation, ft represents the rate of change of the 
spin density in response to an external electric field. 

Intrinsic spin-orbit effects have been shown to lead to 
a non-zero Berry curvature which gives a contribution 
to the anomalous Hall effeclj 2 ^ 2 ^!, while a spin-orbit- 
induced metal-insulator transition has been detected by 
Koga et al£&. Aside from the references mentioned 
above, such intrinsic spin-orbit effects have not been 
taken into account previously, although several theories 
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exist which account for the role intrinsic spin-orbit ef- 
fects play in spin relaxation^. Moreover, a number of 
articles have considered an electric-field-induced rotation 
of the total spin polarization£2i2i, which is distinct from 
the electrical spin generation we discuss in our work. In 
the situations we consider, the total spin polarization is 
initially zero, therefore it cannot undergo a rotation. 

Because both the rate of change of the spin density and 
the electric field are even under time reversal, ferromag- 
netism is not required for spin generation. Nevertheless, 
because the electric field changes sign under spatial in- 
version while the rate of change of the spin density does 
not, spin generation occurs only in crystals with broken 
inversion symmetry. Further details of these symmetry 
arguments will be given in the last section. 

The outline of this paper is as follows. In section II we 
will present the formalism underlying our discussion of 
intrinsic spin generation. The mechanism we outline ap- 
plies to a wide class of systems, and to illustrate this we 
discuss, in section III, two-dimensional heterostructures 
described by the Rashba model, followed by a model of 
the conduction bands of unstrained bulk wurtzite struc- 
tures and strained bulk n-GaAs. 



II. FORMALISM 

In the physical picture we consider, the spin-orbit in- 
teraction has been taken into account in the band struc- 
ture. We adopt a Boltzmann-wave-packet approach, in 
which the carriers in a band labeled by the index n are 
described by wave-packets following a Boltzmann phase- 
space distribution /„(r c ,k, t). The construction of a 
wave-packet representing a charge and spin carrier in 
band n, which has real and A;-space coordinates (r c ,k), 
has been thoroughly treated by Sundaram and Niu2£ and 
will not be considered at length here. The semiclassical 
equations of motion for (r c , k), in the presence of a con- 
stant uniform electric field E, are^: 



Boltzmann equation, 



Kk r = — eE 



de r , 



Hr c = — h eE x (l n , 

ok c 



(1) 



where fl n = 2Im(^ t |^ L ) represents the Berry, or geo- 
metrical curvature^. 

The theory presented in this article is a theory of spin 
accumulation in the strong spin-orbit coupling limit, im- 
plying that the splitting of the bands due to the spin-orbit 
interaction exceeds their broadening due to disorder at 
all wave- vectors. As a result, the bands are well defined 
although, due to the presence of the spin-orbit interac- 
tion, they are not pure spin-up and spin-down bands. 
The distribution function corresponding to each band n, 
f n , can drift according to the semiclassical equations of 
motion, and can also change due to collisions. The time 
evolution of the distribution function is governed by the 
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The right hand side, ^f| co ;/, is the collision term which 
may be modeled by a relaxation time approximation 
when scattering is weak enough, or it may be expressed 
in terms of collision integrals. For electrons in equilib- 
rium, the solution of the Boltzmann equation is f n °^ , the 
Fermi-Dirac distribution, while in a general nonequilib- 
rium situation f n can be written as + Sf n . In the 
presence of a constant uniform electric field E the shift 
5f n is given, in the relaxation time approximation, by 
the well-known result: 



fifn — fn f n 



(0) 



eTpEi ■ v, 



(0) 
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and t p is 



where v„ = jr , e n is the energy of band n 
the momentum relaxation time. We shall refer to terms 
which depend only on the equilibrium value of the distri- 
bution function, f n °\ as intrinsic, as opposed to extrin- 
sic terms, depending on the nonequilibrium distribution, 
and therefore on scattering. 

The study of spin generation necessarily relies on the 
spin equation of continuity. For the case we consider, we 
have shown in a previous publication^ that this equation 
takes the form: 



The terms on the LHS represent the spin density and cur- 
rent in band n, while the last term in the equation takes 
into account collisions. The abbreviation (s) n stands for 
the expectation value in band n of the spin operator cor- 
responding to any one component of the spin. We spe- 
cialize henceforth in homogeneous systems, in which the 
divergence of the spin current in J3J will be zero and the 
equation of continuity may be written as: 



dSn 

dt 
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where r p is the relaxation time. The first term on the 
RHS, which accounts for spin generation in the absence 
of collisions, is the focus of this paper. This term, which 
we shall call the torque density, exists due to the fact 
that spin is in general not conserved and therefore the 
average spin of a wave-packet is not constant in time. As 
discussed in Culcer et alm^, this torque density is defined 
as: 



T n (r,t) 



d 3 r 



d 3 fc/„(r c ,k, t)<fJ(f-r)>, 



(6) 



in which f is understood as i [H, s], H is the Hamilto- 
nian, f is the quantum-mechanical position operator and 
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() n stands for the expectation value in a wave-packet con- 
structed starting from the eigenfunction corresponding to 
band n. Throughout this paper we assume that products 
of non-commuting operators have been symmetrized. We 
note that the equation of continuity @ is derived directly 
from the first-principles definitions of S n , 3 n and %^>. 
In homogeneous systems the torque density simplifies to 

T n = Jd 3 kf n {r) ni (7) 

which will be referred to as the spin generation term. 
The fact that we are considering homogeneous systems 
also implies that we may regard the wave-packets as be- 
ing wide in real space, thus sharp in fc-space, and evaluate 
the expectation value (f }„ using Bloch wavcfunctions. It 
should be pointed out that, although we arrive at our 
results semiclassically, one does not require a local de- 
scription to obtain them, and they can be found using, 
for example, a Kubo formula approach. The connection 
to this latter approach will be detailed in what follows. 
Moreover, our theory is not restricted to the generation of 
spin, since (s) n may represent the wave-packet expecta- 
tion value of any component of any other non-conserved 
observable. 

We will be concerned with a system in which only a 
constant uniform electric field is present. Making a con- 
venient choice of gauge, this electric field can be included 
in the Hamiltonian through the electromagnetic vector 
potential A(r, t) only. This results in a nonadiabatic 
mixing of the bands, so that the Bloch wavefunctions 
\u n ) have the following form: 

■ > I / \ \- (<f>m\ihj- t \(f> n ) 

\Un) = \4>n) - 2^ Z I0™>' ( 8 ) 

where the are the unperturbed Bloch eigenstates. 

The only time dependence comes from the fact that 
k drifts under the action of the electric field, as in 
Therefore it is legitimate to make the replacement 
^ = — ^ ■ ^ in ©. In this way, the wave functions 
\u n ) depend on the electric field through a reactive term, 
in other words the field induces a change in the wave 
functions at each k. The {|w n )} form a complete set. 
Moreover, the expectation value (s) n = (u n \s\u n ) is a 
function of k only, and its time dependence arises im- 
plicitly through its dependence on the wave vector. We 
have thus included the effect of the electric field in mixing 
the bands but neglected inter-band scattering. 

In the limit of wide wave-packets, it is straightforward 
to prove, starting from Eq. JSJ), that (f}„ evaluated in 
the {\u n }} basis is equal to = k • where (s) is 
evaluated in the unperturbed basis and k = — ^ 

from Q). The former approach is equivalent to using 
the Kubo formula to find the response of f to an electric 
field. Following this line of thought, we find that the 
spin generation term is always at least first order in the 



electric field. Then, to first order in the electric field, we 
may replace /„ by its equilibrium value f n °\ This spin 
generation term is then purely intrinsic, as defined at the 
beginning of this section. The final form of this term is: 

Our theory thus shows that there exists a spin generation 
rate which can be interpreted as due to a displacement 
of the wave-packet in fc-space, but also as the expecta- 
tion value of the operator f in a state which is not an 
eigenstate of the crystal Hamiltonian. Experiment*^ has 
been attempting to measure this rate of generation, but 
for the time being success has been achieved only in the 
weak spin-orbit coupling limils^i. 

In the steady state the equation of continuity becomes 
simply: 

I rf a fc (/ "~ /n % )„=T» (10) 

J T p 

Remembering that J d 3 kf n °\s) n — and assuming a 
momentum relaxation time independent of wave-vector, 
the LHS is simply — . We can then rewrite the above 
equation to express the steady state spin density as: 

S n = T n T p (11) 

The characteristic time governing the relaxation of the 
spin density distribution is the momentum relaxation 
time. This is to be expected, since the theory describes 
non-degenerate, well defined bands. 

We note that the torque term must be present even 
in a clean system if the Hamiltonian contains spin-non- 
conserving terms. In the presence of scattering mecha- 
nisms, the intrinsic spin generation term is balanced by 
the extrinsic spin relaxation term so that a net spin polar- 
ization can be reached in the steady state. In addition, in 
the systems we consider, we assume scattering is strong 
enough to keep the distribution function near equilibrium 
and the scattering time small, but not strong enough to 
make inter-band mixing important. 

III. EXAMPLES 

In order to clarify the significance of the examples pre- 
sented below, we start with some comments regarding 
the spin-orbit interaction and asymmetry. The terms in 
the spin-orbit interaction which are odd in fc rely upon 
the inversion asymmetry of the system under study. This 
asymmetry can be of two kinds, depending on the dimen- 
sionality of the system. In three dimensions, the inver- 
sion asymmetry is a property of the underlying material, 
and is referred to as bulk inversion asymmetry (BIA). 
In two dimensions, an asymmetric confinement potential 
can provide an additional source of inversion asymmetry, 
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known as structure inversion asymmetry (SIA). More- 
over, the application of strain along a particular direc- 
tion further reduces the symmetry of the structure, with 
important consequences which will be examined below. 



A. Rashba-type spin-orbit interaction 

Our theory allows us to treat bulk semiconductors and 
quantum wells on the same footing. We begin with a 
study of the Rashba Hamiltonian. This Hamiltonian de- 
scribes, based on symmetry arguments, the SIA of quan- 
tum well or heterojunction based two dimensional elec- 
tron systems and is usually the dominant source of spin- 
orbit coupling in these systems. The effective Hamilto- 
nian has the form: 



H = — 1- a(a x k) • z, 

2m 



in which a is the spin orbit constant and a is the vector of 
Pauli spin matrices, a is usually taken to have the form 
a>A&E z , where is a material specific parameter while 
E z is the component of the electric field in the z direction. 
This electric field is in general a function of position in the 
quantum well, and there exists in principle an additional 
contribution if the interface on one side of the quantum 
well is different from the interface on the other side. This 
effective Hamiltonian therefore provides only an approxi- 
mate description of the Rashba effect. The magnitude of 
the Rashba interaction can be tuned by an external gate 
voltage by an amount which has been shown to be as 
much as 50%2». It is customary to view the term multi- 
plying the spin as a momentum-dependent effective mag- 
netic field 15 in which the spin precesses. In the absence 
of an external magnetic field, the bands in the Rashba 
model are degenerate at k = 0, and each band contains 
the same number of spin up and spin down carriers^ 5 .. 
The Hamiltonian has eigenvalues e± = Ar^- ± ak, which 
will be labeled by + and - respectively. The Berry cur- 
vatures are Cl± = T^limg^o ^^i/; . The spin 
generation term takes the following form: 



(r) ± = T^(kxE) 



z. 



(13) 



with (f ) defined in Eq. Q . Interestingly, the spin torque 
does not depend on the spin-orbit constant. However, 
the total torque term, summed over the two bands, de- 
pends on the difference in Fermi wave vectors, which is 
proportional to the spin-orbit constant. We find that 



_ earn _ 
T =4^ EXZ ' 



(14) 



which vanishes in the limit in which a — > 0, is in agree- 
ment with Magarill 24 >. This result does not depend on the 
number density and has a universal form, but it should 
be noted that its validity is not universal, rather it is re- 
stricted to systems in which disorder is weak as discussed 
at the end of the previous section. 



Using symmetry arguments, we find that the spin-orbit 
coupling in the conduction band of bulk wurtzite struc- 
tures is also described by a Rashba-type Hamiltonian, 
with a spin-orbit constant defined analogously. This con- 
clusion is supported by group theory arguments^ 6 . The 
only terms linear in k (and therefore dominant except 
at high densities) which are allowed by symmetry are 
P{a x k y — a y k x ), and the Hamiltonian is: 



h 2 k 2 

H = — (- (3{a x k y - (jykx) 

2m 



(15) 



with eigenvalues (labeled as before) 4j-^- ± (3k. The spin 
generation term has a form very similar to 113|) : 



(r)± = ^^(kx x E) • z 



(16) 



(12) 

In the above, k± = (k x , k y ,0). The total torque term is: 



47T 2 ^ 2 V ' 



(17) 



which again vanishes as (3 — > and as the number density 
n -> 0. 



B. Cubic Dresselhaus spin-orbit interaction 

Finally, we turn our attention to the conduction band 
of zincblende semiconductors, which has been the focus 
of experiment recently^. In order to be close to exper- 
iment, we consider a degenerate electron gas in an n- 
doped In^Gai-^As heterostructure grown on GaAs, with 
x=0.07, with a strain of 0.46% directed along (001), as 
given in the recent paper of Kato et al^ and references 
therein. The Fermi surface is only slightly displaced from 
equilibrium, as revealed by the drift velocities and by the 
mobility measurements in the supplementary table. We 
will not, however, attempt to simulate experiment, as 
many aspects seem to remain incompletely understood. 
For example, it is known that the strain tensor acquires 
off-diagonal shear terms^I, but experiment does not so 
far provide an unambiguous way of determining their 
rol o 19 i 38 . Therefore, we will consider an idealized situ- 
ation in which the x and y lattice constants are equal to 
their substrate values, while the lattice constant in the 
z direction expands according to the elastic equations. 
Moreover, there is no convincing pattern in the variation 
of the BIA with increasing strain in the experiment. In 
fact, the only sizable increase is observed at the largest 
value of the strain, which is 0.46%. Thus, it is not clear 
whether the linear or cubic term in the BIA is dominant 
in the samples investigated. 

The symmetry of zincblende does not allow terms lin- 
ear in k in the conduction band in the bulk. As a re- 
sult, when strain is applied these linear-in-A; terms will 
be first order in the strain. These terms will be impor- 
tant at small wave vectors, but we will concentrate on 
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situations in which the number density makes the Fermi 
wave-vector kp high enough that the cubic Dresselhaus 
term is dominant. We therefore neglect the effect of the 
terms linear in k in this calculation, although we take 
into account the effect of strain on the effective masses. 
We will take into account only the spin-orbit terms cubic 
in k which are present in the unstrained lattice 3 ^, namely 
Xo x [k x {ky — fef ) +c.p.], where A is the spin-orbit constant 
and c.p. stands for cubic permutations. To determine 
the range of validity of this approach, we estimate the 
doping density n at which the k 3 term dominates the 
linear-in-fc BIA and S1A terms. Based on the data in 
Kato et ali^, we estimate that the term cubic in k will 
dominate if n > 2.7 x 10 16 cm~ 3 , which puts the ex- 
periment, in which n = 3 x 10 16 cm" 3 , narrowly in the 
range in which this term is dominant. Our prediction is, 
however, consistent with the bulk BIA findings of Kato 
et al. 

Our theory is valid in the limit of strong spin-orbit 
interactions or weak disorder. For our theory to be valid, 
the following must hold: 



— < A so (k F ), 

Tn 



(18) 



where t p is the momentum relaxation time and A so (fcj?) 
is the spin-orbit splitting at the Fermi wave vector. From 
the mobility, t p is estimated at 0.22ps, yielding ^- = 2.9 
meV. In order for the spin-orbit to overwhelm this, the 
number density must exceed 3.5 x 10 18 m~ 3 , which is 
a more stringent requirement than the requirement that 
the cubic Dresselhaus term exceed the linear one. In 
other words, when the system is in the strong spin-orbit 
limit it is already in the regime where the fc 3 term dom- 
inates. This calculation also shows that the experiment 
of Kato et al. lies in the weak spin-orbit coupling limit, 
outside the validity limit of our theory. The experiment 
is performed in a doped semiconductor in the extrinsic 
regime, where extrinsic (as defined in section II) refers to 
situations in which scattering effects are dominant, that 
is the spin-orbit splitting at kp is broadened beyond the 
point at which the bands overlap. By increasing n one 
increases hp and therefore A so (kp), passing into the in- 
trinsic regime. 

The Hamiltonian for this system is: 



2 1,2 



H 



K 2 k 



2J.2 



2m j 



h 2 k 

2m. 



+ \(T X [k x {k 2 v - k 2 z ) + c. P : 



(19) 



where kj^ has been defined above. It has eigenvalues 



trk- 



2 — - ± AA (labeled as before), with A 



given by J[k 2 {k 2 — k 2 ) 2 + c.p.]. The Berry curvatures 



are given by Sl± = ± 
model the x-component of spin takes the form: 



(fcg-fe;xfc;-fc;)(fcg-fc> 

2A3 



hk x (k 2 k 2 ) 



2A 



In this 



(20) 



with the other components given by cubic permutations 
of this expression. The rc-componcnt of the spin genera- 
tion term is: 



eXE x (k y k z )(k y k z + k z k y k x k y k x k z ) 
2A 3 



(21) 

Again the other components can be found by cubic per- 
mutation. Note that, if strain were absent so that the 
effective mass would be isotropic, the band structure 
and thus the equilibrium distribution would have cubic 
symmetry, making (f ) zero after integration over k. As 
a result, we expect that in the strong spin-orbit cou- 
pling limit the spin generation will be proportional to 
the strain. 

To facilitate comparison with experiment, it is more 
instructive to examine the electric torque response ten- 
sor, which we define through the equation % — xJjEj, 
and x xx = — Xyy have finite values, the other compo- 
nents being zero. In contrast to the Rashba model, the 
diagonal components of the tensor are finite, whereas the 
off-diagonals vanish. To obtain an explicit expression for 
the total spin torque, we must integrate over the Fermi 
surface, which in this case is ellipsoidal. We use the equi- 
librium distribution and the fact that /(e±AA) = f(e) ± 



AA 



a/(e 



ith e = 



trk. 



K K J 9f(e) c/ \ 



de ' ,vilJi - 1 ^ 2m± 1 2m z 

We write the mass ratio as = 1+7, where 7 is a 
small quantity, and we find that for the system under 
study 7 = 0.023 39 . The total result is: 



T ine\m\T 
Xxx = 8^ 2 ( +l) 



1/2 



(22) 



where X is a dimcnsionless angular integral, which con- 
tains the angular part of 121(1. given by: 



J= f del d<\> sin 6>[sin 2 9 sin 2 4> - (1 + 7) cos 2 9} x 4!r ^ , (23) 
Jo Jo 0(0, W 

where the functions / and g are as follows: 

f(9, 4>) = (1 + 7) 2 cos 4 9 sin 2 + (1 + 7) cos 2 6 sin 2 9 sin 4 4> - sin 4 9 cos 4 (f> sin 2 <j) - (1 + 7) cos 2 9 sin 2 9 cos 4 (f> 
g(6,4>) = sin 4 9 cos 2 sin 2 </>+(! + 7) sin 2 9 cos 2 6>(cos 4 4> + sin 4 (f>) + (1 + 7) 2 cos 4 6 - 6(1 + 7) cos 2 9 sin 2 9 cos 2 sin 2 0. 
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Using a Monte Carlo integration method, we find that 
X=-0. 03037. Based on symmetry arguments, I vanishes 
when 7 = 0, i.e. when m z = mj_. This can be seen 
from Eq. (|21|l by switching k y and k z . At small 7, 
< 7 < 0.1, our calculations show that T is linear in 
7 and is given approximately by X = — 1.287. We expect 
that, for n = 10 19 cur 3 , X T XX = -1-15 x 10~ 7 c/m 2 . This 
result is four orders of magnitude larger than that ob- 
served in the current experiments in the weak spin-orbit 
coupling limit, which should offer an incentive for doping 
the samples in order to move into the strong spin-orbit 
limit. Meanwhile, for the density n — 3 x 10 16 cm" 3 
used in experiment we find x xx — —3.83 x 10~ 10 c/m 2 , a 
number that is one order of magnitude higher than that 
observed. We stress again however that the experiment 
was performed in the weak spin-orbit coupling limit. 

IV. SYMMETRY CONSIDERATIONS 

The form of the electric torque response tensor in dif- 
ferent systems can be understood based on symmetry ar- 
guments. In Table I we have listed a number of possible 
spatial transformations which are relevant to our prob- 
lem as well as the behavior of the electric field and spin 
torque under these transformations. The transformations 
considered are not assumed to be symmetry operations of 
the materials under study. I m refers to spatial inversion 
along the m-axis, that is m — > — m, while R m refers to a 
rotation of arbitrary magnitude anticlockwise about the 
m-axis. For simplicity and without loss of generality we 
take the electric field to be directed along x and consider 
the generation of spin along all three Cartesian axes in 
response to this electric field. 

Referring to the first part of the table, first row, it can 
be seen that under spatial inversion in the x-direction the 
electric field changes sign while (f) x remains the same. 
On the other hand, under the same transformation (r) y 
and (f) z change sign, behaving in the same way as the 
electric field. As a result, if spatial inversion along a; is a 

TABLE I: Behavior of the electric field and torque under var- 
ious spatial transformations 



Operation 


E x 




(r)y 




I* 




+ 






h 


+ 




+ 




h 


+ 






+ 


Rx 


+ 


+ 






Ry 






+ 




r: 








+ 


T R n /' 2 
^x^x 




+ 




+(r)y 


r b*/ 2 

1 y J^ x 


+ 






-(r)y 




+ 









symmetry operation of the material the torque term (f) x 
will vanish, but no information can be inferred about 
(f)j, and (f) 2 . An analogous argument holds for inver- 
sion along the y and z axes. Therefore if spatial inver- 



sion along all three axes is a symmetry of the material all 
components of the torque must vanish. This is consistent 
with the observation that (f } = under a three dimen- 
sional inversion. Furthermore, this table helps explain 
the form of the torque in systems with Rashba spin-orbit 
interaction. The Rashba Hamiltonian has I x and I y as 
symmetries while I z is broken. Examining Table I we see 
that in response to E x only (r) y can be nonzero, consis- 
tent with our findings in the previous section. 

Next let us consider the effect of rotations by tt about 
the three Cartesian axes. In the second part of the table, 
first row, we examine a rotation about x. This rotation 
does not affect E x and (f) x , but (f) y and (f) z change 
sign. Therefore, if R x is a symmetry operation (r) y and 
(t) z will be zero in response to E x . By extension, if rota- 
tions by 7r about all Cartesian axes are symmetry opera- 
tions then all the off-diagonal components of the electric 
torque response tensor are zero. For example, in GaAs, 
with or without strain applied along the z-axis, rotations 
by 7r about all three axes are symmetry operations so the 
off diagonal components of \ T are zero. 

Let us also consider the combined effect of rotation 
and inversion, taking GaAs as an example. GaAs does 
not have inversion along any of the three Cartesian axes 
as a symmetry operation. Therefore we cannot infer the 
vanishing or survival of the electric torque response ten- 
sor in GaAs based on considerations of inversion alone. 
With the electric field still along x, we consider an an- 
ticlockwise rotation by \ about the x-axis followed by 
inversion along x, y and z. The rotation does not af- 
fect E x and (t) x . Looking at the third part of the table, 
first row, the electric field changes sign under I X R X while 
(7% does not. Therefore, if 1 X R*J 2 is a symmetry of 
the system, the cc-component of the torque in response 
to E x will vanish. A similar argument holds for the y 
and z axes. Therefore if I m Rm 2 is a symmetry for all 
Cartesian axes all the diagonal components of the elec- 
tric torque response tensor will be zero. This operation 
is a symmetry of bulk GaAs and in order to remove the 
rotational symmetry we consider strain applied along the 

z axis. In this case I X R X ^ 2 and I y Ry^ 2 are not symme- 
tries anymore, allowing x xx an d Xyy t be nonzero, while 
X T ZZ is zero. 



V. CONCLUSIONS 

We have shown that in crystals with inversion asym- 
metry and strong spin-orbit interactions a spin accumu- 
lation will be generated in the presence of an electric field 
due to an intrinsic spin torque term. In the steady state, 
the spin accumulation is given simply by the product of 
this torque term with the momentum relaxation time. 
This term is expected to be observable both in systems 
with Rashba spin-orbit interactions and with Dressclhaus 
interactions. In the latter we have shown that by doping 
so as to bring the system into the intrinsic regime, the 
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effect can be observed using currently available technol- 
ogy. 

Finally, we would like to mention a novel experimental 
method which has been recently proposed. In condensed 
matter it is usually not possible to isolate a single charge 
and spin carrier in order to measure effects occurring on 
the scale of an individual wave-packet. Therefore, to be 
able to verify the existence and properties of the spin gen- 
eration term one may resort to atomic physics. Using a 
cold-atom system to construct an individual wave-packet 
and mimic the spin-orbit interaction, an experiment can 



measure, for example, the size, wave vector and electric 
field dependence of (f) for this wave-packet. This method 
and the physics underlying it are described in detail by 
Dudarev et alm&. 
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